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A Neutron Time-of-Flight (NTOF) spectrometer, based at the Heavy Ion Research Facility in Lanzhou 
(HIRFL) was developed for studies of neutron production of proton induced spallation reactions related to the 
ADS project. After the presentation of comparisons between calculated spallation neutron production double- 
differential cross sections and the available experimental data, a detailed description of the NTOF spectrometer 
is given. Test beam results show that the spectrometer works well and data analysis procedures are established. 
The comparisons of the test beam neutron spectra with those of GEANT4 simulations are presented. 
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I. INTRODUCTION 


The interest in spallation reactions [1, 2] has been increas- 
ing with the development of Accelerator Driven Systems 
(ADS) and the applications of Spallation Neutron Sources 
(SNS) in many fields, such as the transmutation of nuclear 
waste, nuclear energy generation, neutron sources for materi- 
al irradiation, and neutron scattering science. The spallation 
reaction is defined as interactions between a light projectile 
(e.g. proton) with GeV range energy and a heavy target nu- 
cleus, which is split to a large number of hadrons (mostly 
neutrons) or fragments. Many nuclear models, such as the 
intra-nuclear cascade-evaporation (INC/E) [3] model and the 
quantum molecular dynamics (QMD) model [4], have been 
developed to study spallation reactions. These nuclear mod- 
els are coded for thin-target simulations. In order to perform 
calculations for practical applications, which involve complex 
geometries and multiple composite materials, nuclear models 
need to be embedded into a sophisticated transport code. Nu- 
clear models are combined with a Monte Carlo transportation 
code, like MCNP [5], GEANT4 [6, 7], and FLUKA [8, 9], 
nucleon meson transport codes (NMTC) [10] are widely uti- 
lized for designing the facilities of engineering applications 
of the spallation reaction. 

However, the nuclear models embedded into the trans- 
portation codes need to be validated by experimental mea- 
surements. Several laboratories [11-14] worldwide have 
constructed specific experimental setups for measuring the 
double-differential cross sections and the spectra of neutrons 
produced in proton induced spallation reactions. In their ex- 
periments, a time-of-flight technique with organic liquid scin- 
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tillators (or plastic scintillators) or/and a recoil proton mea- 
surement combined with a magnetic spectrometer has been 
utilized for neutron measurements. Those experimental mea- 
surements have made great contributions to the improvement 
of nuclear models. A Neutron Time-of-Flight (NTOF) spec- 
trometer, which is based at the Heavy Ion Research Facility 
in Lanzhou (HIRFL), is designed for further investigation of 
reaction mechanisms and neutron productions of proton in- 
duced spallation reactions related to the ADS project. 

In this work, the prediction ability of three nuclear models 
embedded in the GEANT4 for double-differential cross sec- 
tions of spallation neutron production is examined by com- 
paring the calculated results and the available experimen- 
tal data [15]. The configuration and test beam results of 
the NTOF spectrometer are also presented and discussed in 
detail. 


II. THEORETICAL CALCULATIONS 


GEANT4 is a Monte Carlo transport code developed in 
CERN for simulating the passage of particles through mat- 
ter. It is widely used in particle and nuclear physics, accel- 
erator physics, medical science, astrophysics, and aerospace 
studies. In GEANT4, users have abundant choices of physics 
models to handle the interactions of particles with matter over 
an extended energy range. 

Spallation reactions are usually described in two stages: 
the intra-nuclear cascade and de-excitation. At the first stage, 
the incident particle transfers its kinetic energy to target nu- 
cleons by elastic collisions and a cascade of nucleon-nucleon 
collisions. Some of the particles that obtained enough energy 
to escape from the nucleus are emitted. The residual nucle- 
us before the de-excitation stage can emit particles with the 
low energy, which are called pre-equilibrium particles, be- 
fore the de-excitation stage. The energies of pre-equilibrium 
particles are generally greater than the energies of particles 
emitted during the later stage. In the de-excitation stage, the 
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Fig. 1. (Color online) Comparisons of experimental neutron pro- 
duction double-differential cross sections for 800 MeV proton on a 
Tungsten target with simulated results of different nuclear models at 
30° to 150° detection angles. 


excited nuclei lose their energy through the evaporation of 
neutrons, protons, or light charged particles. If the nuclei do 
not have enough energy to evaporate the particles, it may emit 
gamma rays and decay to a stable state. 

In this work, the BIC [16], INCL [17], and QMD [18] 
models embedded in GEANT4 were used to calculate the 
first stage of the reactions. A statistical model is used for 
the second stage. The dynamical model is linked to the de- 
excitation handler provided by GEANT4. As a result, the 
double differential cross section of neutron production can 
be given at each stage. In Fig. 1, the experimentally mea- 
sured neutron production double-differential cross sections 
from the reaction on the Tungsten target induced by 800 MeV 
protons at detection angles of 30° to 150° are compared with 
the simulated results of BIC, INCL, and QMD models em- 
bedded in the GEANT4. The experimental data were tak- 
en from the EXFOR database [15]. Similar comparisons are 
given in Fig. 2 for 800 MeV to 1600 MeV protons at a 30° de- 
tection angle. Both of the spectra are multiplied by a factor of 
107” (n = 0, 2, 4, 6) from top to bottom. Good agreements 
are achieved among the INCL calculations and experimental 
measurements in the entire angular and energy range exam- 
ined. The BIC and QMD calculations slightly underestimate 
the experimental results in some parts of the neutron spec- 
trum, especially in the high energy tail. These discrepancies 
are related to the neutron detection angle and the incident pro- 
ton energy. 


II. NTOF SPECTROMETER 


The NTOF spectrometer is capable of measuring neutrons 
in wide energy and angular ranges. It consists of a beam pick- 
up detector and 10 individual neutron detection modules. The 
schematic view of the NTOF spectrometer is shown in Fig. 3. 
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Fig. 2. (Color online) Comparisons of experimental neutron pro- 
duction double-differential cross sections for 800 MeV to 1600 MeV 
protons on Tungsten target with simulated results of different nucle- 
ar models at 30° detection angle. 
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Fig. 3. A schematic view of the NTOF spectrometer. 


The beam pick-up detector is composed of a BC404 plastic 
scintillator (5cm x 5cm x 0.2 cm detector with a dual-PMT 
readout at both ends and is located upstream from the tar- 
get. It is used to obtain the time and position of the inci- 
dent particle beam. Each of the neutron detection modules 
is composed of a thin plastic scintillator detector (veto de- 
tector) and organic liquid scintillator detector. Two different 
sizes of cylindrical BC501A organic liquid scintillator detec- 
tors are placed at different angles. The larger BC501A detec- 
tors have a 12.7 cm diameter and are 12.7 cm in length. They 
are used for detecting higher energy neutrons with longer 
neutron flight paths. The smaller BCS501A detectors have 
a size of 5.08cm in diameter and 5.08cm in length and 
they are used to measure low energy neutrons. In front of 
the individual BCSO1A detectors, a BC404 plastic scintilla- 
tor (15cm x 15cm x 0.3 cm) coupled to a 9813KB PMT are 
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Fig. 4. (Color online) Separation of charged particle events from 
non-charged particle (neutron and gamma-ray) ones, using the veto 
counter and a BC501A liquid scintillator. 


mounted as a veto detector to distinguish charged particles 
from non-charged particles (neutrons and gamma-rays). For 
the separation of gamma-rays from neutrons, the pulse shape 
discrimination (PSD) property of the organic liquid scintilla- 
tor detector is utilized. Neutron kinetic energy is calculated 
from the time-of-flight (TOF) spectrum between the organ- 
ic liquid scintillator detector and the beam pick-up detector. 
The neutron flight path from the target to neutron detector is 
about 1.5m. The 400 ns TDC range is used in the experiment 
for measuring above 0.1 MeV neutrons. 

A VME-based data acquisition system (DAQ) for record- 
ing the experimental data on an event by event basis has been 
developed. The DAQ software is based on the CERN ROOT 
framework and runs on a Linux operating system. The main 
functions of the DAQ include data acquisition, online moni- 
toring and offline data analysis with a graphical user interface 
(GUD), which provide convenient operation. 


IV. RESULTS AND DISCUSSION 


A test experiment for the whole system of the NTOF spec- 
trometer was performed by measuring neutrons, gamma rays, 
and charged particles in the energy spectrum, production 
yield, and angular distribution from a Tungsten target with 
an 16O beam at a bombarding energy of 400 MeV/u. 

In the experiment, the combination of the light output spec- 
tra of the veto detector and the BCSOIA scintillator detec- 
tor was utilized for the separation of charged particle events 
from non-charged particle (neutron and gamma-ray) ones, as 
shown in Fig. 4. The anode signal of the BCS01A scintillator 
detector was divided into three pulses. Two of them were fed 
into charge-to-digital converters (QDCs), which had different 
gate widths (total and slow). QDCs with total and slow gates 
were used to eliminate gamma-ray events with PSD by the t- 
wo gate integration method. A typical result of PSD is shown 
in Fig. 5. The TOF spectra were measured between the beam 
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Fig. 5. (Color online) A typical result of the pulse shape discrimina- 
tion (PSD) in a BCSO1A liquid scintillator. 
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Fig. 6. (Color online) Time-of-flight (TOF) spectrum for neutrons 
and gamma rays. 
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Fig. 7. (Color online) Comparisons of experimental neutron produc- 
tion yield with GEANT4 results for 400 MeV/u '°O bombarded on 
Tungsten at detection angles of 10° to 60°. 
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pick-up detector and the BCSO1A detector. In Fig. 6, a typ- 
ical TOF spectrum, excluding the charged particle events, is 
given. The prompt gamma peak of the TOF spectrum was 
used as a time reference of the spectrum. The energy calibra- 
tion of the organic liquid scintillator detectors was accurately 
determined by comparing the experimental light output of s- 
tandard gamma sources with a GEANT4 simulated one [19]. 
In Fig. 7, the comparisons of the experimental neutron pro- 
duction yield at detection angles ranging from 10° to 60° us- 
ing lager BCSO1A detectors with GEANT4 results are shown. 
The neutron production yield is converted from TOF spectra 
with normalizing by unit solid angle and incident ion num- 
bers. The experimental neutron spectra shape are well repro- 
duced by the simulations. The experimental data have been 
normalized to the simulated one at 10°. The experimental re- 
sults show that the whole system of the NTOF spectrometer 
works well, and the data analysis procedure is established. 


V. CONCLUSION 


For nuclear engineering facilities, sophisticated simula- 
tion tools based on nuclear reaction models are required. In 
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